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Abstract: Simulating the k(800) on the lattice is a challenging task that starts to 
become feasible due to the rapid progress in recent-years lattice QCD calculations. As 
the resonance is broad, special attention to finite-volume effects has to be paid, because no 
sharp resonance signal as from avoided level crossing can be expected. In the present article, 
we investigate the finite volume effects in the framework of unitarized chiral perturbation 
theory using next-to- leading order terms. After a fit to meson- meson partial wave data, 
lattice levels for ttK scattering are predicted. In addition, levels are shown for the quantum 
numbers in which the cr(600), / (980), a (980), 0(1020), If* (892), and p(770) appear, as 
well as the repulsive channels. Methods to extract the k(800) signal from the lattice 
spectrum are presented. Using pseudo-data, we estimate the precision that lattice data 
should have to allow for a clear-cut extraction of this resonance. To put the results into 
context, in particular the required high precision on the lattice data, the <t(600), the P- 
wave resonances if* (892) and p(770), and the repulsive 7rK, tttt phases are analyzed as 
well. 
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1 Introduction 



Scalar resonances attract increasing interest as, with the rapid progress of lattice QCD 
simulations [1], the excited resonance spectrum starts to come into reach [2-17], including 
first results on the k(800) [17]. The strangeness S = — 1, k(800) resonance is of particular 
interest because it does not mix with scalar glueballs and is thus expected to reveal itself 
with a clearer signal than, e.g., the broad <t(600) in the scalar-isoscalar sector. On the 
other hand, the k presents a similar challenge, both for experimental and lattice studies, 
as this resonance is similarly broad. For wide resonances, the avoided level crossing of 
the lattice levels is washed out 1 and a clear-cut extraction of the resonance pole and its 
properties requires additional effort. Such complications from finite volume effects can be 
addressed in the framework of field theory with periodic boundary conditions, as pioneered 
by Liischer [19-21]. On one hand, this framework allows to predict lattice levels using 
hadronic approaches originally set up for the infinite volume limit. On the other hand, 
phase shifts, and in consequence resonance properties, can be extracted from existing lattice 
data. In the absence of such data, as it is the case for the k(800), the formalism can still 
be used to study the propagation of the unavoidable uncertainties from lattice data to 
the extracted resonance properties. Predicting lattice levels in the meson-meson sector, 
using the chiral expansion up to next-to-leading order (NLO) terms, and determining the 
needed precision of the lattice spectrum to make sensible statements about the k(800), is 
the purpose of this study. 

For the prediction of lattice levels and the generation of pseudo-data, the chiral uni- 
tary scheme of ref. [22] is employed that includes terms up to NLO in the chiral expansion. 
Chiral Perturbation Theory (CHPT) [23, 24] has been applied to the problem of ttK scat- 
tering since a long time, calculating threshold parameters [25], one-loop corrections [26], 
and extensions to higher energies by including resonances [27]. The unitarization of the am- 
plitude can induce strong non-linearities leading to non-perturbative features such as pole 
formation and has been explored, e.g., in ref. [28] where the <r(600), /o(980), and ao(980) 
resonances appear dynamically generated from the unitarization of the lowest order (LO) 
chiral meson-meson interaction in a Bethe-Salpeter type of equation. This approach has 
been further developed and refined by, e.g., taking resonances into account to extend the 
applicability to higher energies [29-32]. An alternative scheme (among other approaches 
such as meson exchange models [33, 34]) is provided by the inverse amplitude method — the 
approach employed here [22, 35-37]. It was extended to provide the full one-loop calculation 
for the chiral transition potential in SU(3) [38] and to study the quark-mass dependence of 
resonances [39-41]. Dispersive analyses of meson- meson scattering as formulated, e.g., in 
Roy equations provide a related model-independent approach, often stabilized by the use 
of chiral perturbation theory [42-47]. 

To learn about the k(800) resonance on the lattice, the mentioned approaches, all 
of them set up for the infinite volume limit, can be extended to predict lattice levels 
using Liischer's formalism. The idea of using effective field theory to study the lattice 
spectrum and excited states has been been developed for the A(1232) resonance and other 

1 This phenomenon was first studied in Ref. [18] in the context of extracting the A-resonance properties. 
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meson-baryon systems [18, 48-51] as well as for the meson sector [52-54]. A special role 
plays here the generalization of the Liischer formulation to coupled channels [50, 53, 54]: 
e.g., with the quantum numbers of the k(800), one has the 7r.fr and the r]K channels, 
the latter not being too far away from the k(800) pole, thus potentially influencing its 
properties. In such situations, phases and resonances can still be extracted using a K- 
matrix approach [53] together with the Liischer Zqq function or directly the relativistic 
meson-meson propagator [54]. Compared to the one-channel case, for which the Liischer 
formula is valid [19, 20], the phase extraction in the two-channel case requires in general 
much more precise lattice data as recently shown in ref. [54], in particular if a particle 
threshold coincides with a resonance, as, e.g., the /o(980) (see also ref. [55] for an example 
in the charm sector). Concerning the extraction of higher partial waves (see, e.g., refs. [48, 
56]) in the presence of multiple channels, another possible analysis tool for lattice data is 
given by dynamical coupled-channel approaches [57] that were recently shown to be easily 
adaptable to the finite volume problem [58]. 

In the present study, lattice levels will be constructed and analyzed following the 
discretization procedure of ref. [54]. The relativistic propagator G used in that unitary 
chiral (UCHPT) approach leads to a discretized propagator G that provides a formalism 
equivalent to the Liischer approach, plus relativistic effects of order e~ * L , the size of which 
have been discussed in ref. [54]. Apart from these particular e~ M ^ L effects, there are of 
course, effects of the same order from, e.g., omitted loops in the t-channel or polarization 
effects of the pion [21]. Note that the propagator G used in ref. [22] is just the same as 
that of ref. [54]. The present work provides an extension of methods developed in ref. [54] 
to the irK sector and goes beyond it in several aspects: first, the use of NLO terms in the 
chiral expansion implies higher powers of the scattering energy E in the potential. In fits to 
pseudo-data generated from this interaction, the non-zero higher order terms will present 
a new challenge for the resonance extraction and new tools are required. The key is here 
that the well-known LO chiral interaction can be used as explicit input in the expansion 
of the general fit potential, and the unknown higher order effects can be cast in terms of 
a systematic expansion in powers of s = E 2 . As will be shown, the chiral input helps to 
stabilize the extraction process considerably. Also, we formulate a criterion of when to cut 
this expansion in s in the resonance extraction process, guided by statistical arguments. 

It should be stressed that the proposed reconstruction method of the infinite volume 
limit is model-independent in the sense that it does not make assumptions on the existence 
and form of the resonance propagator as, e.g., in refs. [4, 9, 15, 17], sometimes in form of 
an effective range formula [4, 9]. In particular, the low-lying scalar resonances are known 
to be very broad [59-64] and can by no means be approximated by Breit-Wigner or related 
functions. 

In section 2 we briefly introduce the inverse amplitude method of ref. [22] and show 
the result of our fit to the meson-meson partial waves with special weight on the (rjK, irK) 
system in which the k(800) and the if* (892) reside. The lattice levels, predicted from this 
solution, are shown in section 2.1. Having set up an amplitude that describes well the ttK 
scattering data, with low-energy constants well constrained from various different channels, 
the main issue of the paper can be addressed in section 3: depending on the precision of 
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the lattice data, how well can the properties of the k(800) resonance be determined? After 
a discussion of a suitable expansion potential in section 3.1 — see the above remarks — this 
question is addressed in section 3.2, with pseudo-data derived from the global fit obtained 
before. As the pole position in this fit is rather at the higher end of the accepted pole 
positions [59], we will analyze in section 3.3 the opposite case with a rather light k(800) as 
obtained, e.g., in the data analysis using Roy-Steiner equation of ref. [45]. As the required 
precision of the lattice data, to determine the properties of the k(800), is quite high, one 
may also wonder to which extent one can at least make qualitative statements about its 
properties, using only one data point from the lowest level. This question is addressed in 
section 3.4. To put the results of the k(800) into context and compare to the required 
precision for other quantum numbers, section 4 is dedicated to the analysis of the cr(600) 
resonance. In particular, we check whether the required precision as found in ref. [54] is 
modified in the present formalism in which NLO terms play an important role. In the same 
section, we compare the results also to an analysis of the if* (892) and p(770) resonances 
as well as to the repulsive ttK and irir phase shifts. The residues of all analyzed resonances 
are collected in section 4.4. We briefly summarize our conclusions in section 5. 

2 The pion-kaon amplitude in the infinite volume limit 
with next-to-leading order terms 

Organizing the LO and NLO contact terms in the inverse amplitude method as derived 
in ref. [22], and projecting the interaction to S- and P- waves, a large body of partial 
wave data can be described after a fit of the low-energy constants L\ to Lq + L%. This 
successful approach has been refined over the years, see, e.g., ref. [38] including the full 
one-loop calculation and later also two loop corrections [41] . In this study we use the earlier 
treatment of ref. [22] that is sufficient for the present purposes. Note also that the loops 
in the t-channel will only provide contributions to the potential that are exponentially 
suppressed [e~ Mlt L effects) and thus expected to be of minor importance in the study of 
finite volume effects (for sufficiently large M n L). 

The meson-meson amplitude in the formulation of ref. [22] is 

T = V 2 (V 2 -V4-V2G V2T 1 V 2 (2.1) 

with partial-wave projected transition potentials V 2 and V4 from the LO and NLO chiral 
Lagrangian. The meson-meson loop function G is factorized from the potentials V (on- 
shell reduction), so that V and G in eq. (2.1) are simply matrices in channel space. The 
expressions for the V can be found in ref. [22] 2 or easily obtained from there using crossing 
symmetry, see, e.g., ref. [32]. After partial wave projection, the various particle channels are 
combined to the possible quantum numbers in meson-meson scattering. The normalization 
of T in eq. (2.1) is such that tan 25 = — ~ReT/(ImT + 4irE/q c . m ,) with 5 the phase shift and 
q c . m . the momentum of the particles in the center-of-mass frame. See refs. [22, 28] for the 
two-channel case. Fitting the low-energy constants appearing in V4, a good description of 

2 There is a typing error in the denominator of eq. (B16) of ref. [22] that should be 3F%F^. 
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Figure 1. Solution for the meson-meson interaction [(red) solid lines]. For comparison, the results 
of ref. [22] are shown with the dashed lines. Data: for 6®: (pink) diamonds [67], (blue) squares [68], 
(black) circles: average as defined in ref. [29]; for Sy 2 : circles [69], triangles up: average as defined 
in ref. [29]; for S^ 2 : circles [70], triangles up [71]; for S 2 - crosses [72], squares [73], triangles up [33]; 
for Sy 2 : triangles up [74], circles [71]; for 5{: squares [75], circles [76]. 

the data can be obtained. We have reproduced the result of ref. [22] as shown in figure 1 by 
the dashed lines. For the calculation, we use the meson decay constants F w = 92.4 MeV, 
and, from ref. [24], Fk = 1-22 F n and F„ = 1.3 F n to be able to compare with the results 
of ref. [22]. As we have tested, changes from using more modern values of Fk = (1.192 ± 
0.007) F w [65] and F„ = 118.4 ± 8.0 MeV [66] can be absorbed in the low-energy constants 
without modifying the result significantly 
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Figure 2. Further results (red (solid) lines; data not included in x 2 ) for the {ILS) = (0, 0, 0) sector 
[/o(980) meson] (left and center panel) and for the wrj distribution [ao(980) meson] (right panel). 
Data: left and center: diamonds [67], crosses (X) [77], circles [78], triangles [79], squares [76]; right: 
ref. [80], background from there incoherently subtracted, see ref. [29]. 

The agreement of the solution of ref. [22] (dashed lines in fig. 1) with the phase shift 
data is remarkable, except for the combination of isospin, angular momentum, strangeness 
(J, L, S) = (1/2, 0, —1), i.e. precisely the quantum numbers of the ac(800) we are interested 
in. Thus, we have performed a refit to the data shown in figure 1 (i.e. from the respective 
thresholds up to E = 1.2 GeV), with details given in the following: the k(800) channel is 
included in the x 2 with a special weight of 7. The pole position zq of the c(600) in the 
(/, L, S) = (0, 0, 0) channel has been determined very precisely and model-independently in 
ref. [43], zq = (441 ± i 272) MeV. Thus, we simply consider this value as an additional data 
point that is included in the x 2 ■, with an artificial error of \Azq\ = 30 MeV. While the data 
in the k channel are included with an increased weight, the x 2 ' s from the repulsive 1 = 2 
irir and the 1 = 3/2 nK channels are included with a weight of 0.2, while the data from the 
p channel enter only with a weight of 0.1 (note the data in that channel are quite precise). 
The value of the low-energy constant L3 is particularly sensitive to the p(770)-channel 
and has been fixed in the fit to reproduce the bulk features of this resonance. With these 
trade-offs to achieve a better data description in the k(800) channel, the fit was performed, 
paying special attention to avoid spurious poles in the solution — the inverse amplitude 
method in the form used here is, to some extent, prone to generating narrow poles that 
are not physical. An example, from the solution of ref. [22], is visible in the dashed line 
in figure 1, for the 5y 2 phase at higher energies. In the present fit, such solutions were 
discarded. 

The result of the fit is shown by the (red) thick solid lines in figure 1. In figure 2 
the performance of the present fit is compared to additionally available data that were, 
however, not included in the x 2 ■ Part of the data enter the x 2 with a weight different from 
one, and there is not much sense in calculating the reduced x 2 ■, but the quality of the fit is 
obviously good, in particular in the k(800) channel. The trade-off paid for this agreement, 
visible in the /?(770) and the repulsive irir, ttK channels, is relatively small. 

The values for the low-energy constants and q max (the cut-off in the loop function 
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Li 

n C7Q+0.017 

u.o/o_ 028 


L 2 

627+ ' 028 
U - DZ ' -0.014 


L 3 

-3.5 [fixed] 


La 

-0 7in+ - 022 

u -' iu -0.026 


9 OQ7+0.048 
z - yo ' -0.094 


Lq + L$ 

1 qor+0.026 
1-JOD_0.050 


L 7 

+0.106 
u - '^-0.074 


gmax [MeV] 

981 [fixed] 



Table 1. Fitted values for the Li [xlO 3 ] and q max [MeV] with their uncertainties. As mentioned 
in the text, L3 and g max have been fixed in the determination of the error. 
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pu 


434+i 261 


-31-H9 (KK) 


-30+i 86 (vrvr) 
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[411+i 233] 


[0.2-i 0.2 {KK)} 


[-29+^42 (vrvr)] 











/o(980) 


pu 


1003+i 15 


16— z 79 (KK) 


-12+i4(7T7r) 
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-1 


k(800) 


pu 


815+i226 


-36+i 39 (r]K) 


-30+i 57 (ttK) 










uu 


[804+i 103] 


[-0.4+* 15 (rjK)] 


[-0.1+t 34 (ttK)] 


1 








a (980) 


pu 


1019-i4 


-10-tl07 (KK) 


21-131 (7T7/) 





1 





0(1020) 


p 


976+iO 


-2+iO(KK) 




1/2 


1 


-1 


IT* (892) 


pu 


889+^25 


-10+i0.1 (rjK) 


14+i4(vrK) 










uu 


[915+« 22] 


[-6.4-H0.1 (rjK)] 


[12+i 3 (ttK)] 


1 


1 





p(770) 


pu 


755+z 95 


-U+i2(KK) 


33+il7(7T7r) 










uu 


[791+i 80] 


[-7+il(KK)} 


[27+i 13 (vrvr)] 



Table 2. Pole positions zo [MeV] and residues a_i [M ff ] in the global fit to meson-meson observablcs. 
I, L, S denote isospin, angular momentum, and strangeness, respectively. For the sheets, e.g. pu 
means physical (unphysical) sheet for the channel with higher (lower) threshold, in order. Values 
in brackets refer to poles that should not be compared to PDG values, because they are on hidden 
sheets. For the residues, e.g. (irK) means the residue of T(E) of eq. (2.1), of the irK — > irK 
transition. Note that for the poles in the lower half plane, ci-i(zq) — a*_i(za). 

G [22]) are shown in table 1. These numbers are similar to those of ref. [22] and also of 
the same order as they arise in CHPT, although there is no need for them to be precisely 
equal to a standard CHPT calculation due to a different fitted data range and conceptual 
differences to a full 0(p 4 ) calculation. In addition, the nondinear parameter errors of 
the fit were determined and are also quoted in table 1. These errors are only indicative, 
because there are data in the fit that were included in the x 2 with weights different from 
one; the parameter errors are probably considerably larger than those quoted in table 1. 
Also, in the determination of the error, L% and. {/max 

needed to be fixed. Otherwise, too 
large parameter correlation would have inhibited the determination of the error using the 
usual Xbcst + 1-criterion, where Xbest ^ s ^ ne ^ a ^ ^ s iidnimum. Apart from this, the errors 
indicate that the low-energy constants can be well constrained from the data. 

The resulting pole positions and residues are shown in table 2. If there is more than 
one entry for a resonance, the one without brackets corresponds to the standard sheet that 
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is considered for the resonance. For example, the cr(600) resonance is on the physical sheet 
with respect to the heavier KK channel, but on the unphysical one with respect to the 
lighter irir channel, i.e. on the pu sheet. The uu sheet with the other pole of the a is 
far away from the physical axis at the a energies. For a general discussion of the analytic 
structure in hadron-hadron scattering, see, e.g., ref. [81]. The ao(980) is special because 
its pole is on the (pu) sheet, but at an energy above the KK threshold. This means the 
closest point of the physical axis to this pole is the KK threshold itself, which indeed 
shows a pronounced cusp in the irr] invariant mass spectrum as shown in figure 2. In the 
solutions of refs. [22, 28] the situation of the ao(980) pole is qualitatively the same. The 
role of threshold cusps and the possibility to disentangle their structure from lattice levels 
has been discussed in ref. [54]. 

The pole position of the o"(600) shown in table 2 is close to the value extracted in the 
dispersive analysis of ref. [43] which is not a surprise as that pole position is included in 
the x 2 - The pole position of the K*(892) is very close to the value quoted in PDG [59]. 
This is important as the main objective of this study is the (r]K, ttK) system. 

We would finally like to stress that in the presently considered UCHPT approach the 
considered channel space is limited to the two-body channels of scalar mesons. Multi- 
meson states or channels with (axial) vector particles are omitted, mostly because there 
is no phenomenological requirement for their inclusion at the energies of interest. The 
finite volume structure of multi-meson states is considerably more complicated than the 
two-body case considered here, and beyond the scope of this study. 

2.1 Predicted lattice levels in the meson-meson sector 

The discretization of the loop integral G in eq. (2.1) generates the finite volume effects. In 
ref. [54], the discretized meson- meson loop function for S- waves has been derived with the 
result 



G j = G s,j + (Gj - G s , 



3> ' 



| Qmax | 

q - L n, n X , J M ) - ^^^^ ^ _ ( Ujtl ( q ) + Uj ^[ He ] ^ 



with ujjj = y Mf + |q| , Mj the masses of the two mesons in channel j, and F is a form 
factor introduced in ref. [54] to avoid artefacts (e~ effects) from the sharp cut-off at 
qmax- The function Gj is independent of this form factor, and Gsj — Gsj can be rewritten 
in terms of the Zqq function [54] up to relativistic effects that are of order e~ * . As has 
been studied in ref. [54], the relativistic effects provided by eq. (2.2) can become significant. 
Up to the same e~ Mn L effects the present discretization procedure is equivalent to the K- 
matrix formulation derived in ref. [53]. Note also that in the extraction of the scattering 
phase, in the one-channel case only the combination Gj — Gj enters that is independent of 
the used form factor and cut-off at q max [54] . 
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The lattice levels are obtained from the singularities of eq. (2.1) with the discretized 

Gj, 

&et\V 2 -V A -V 2 GV 2 [ =0, G = dmg{G l ,G 2 ) (2.3) 

where G is a diagonal matrix in channel space that has two dimensions for the cases 
(I, L, S) = (0,0,0), (1,1,0) [KK, vrvr channels], (I, L, S) = (1/2,0,-1), (1/2,1,-1) [rjK, 
ttK channels], and for (I,L,S) = (1,0,0) [KK, ttt] channels]. Only one channel is present 
for the cases (I,L,S) = (3/2,0,-1) [ttK channel], (I,L,S) = (2,0,0) [kit channel], and 
(I, L, S) = (0, 1, 0) [KK channel], in which case eq. (2.3) is one-dimensional with G = Gj. 

For the P-wave states, one has to complement eq. (2.2) with the corresponding angular 
momentum structure for the vertices, as spherical symmetry is broken on the lattice. In the 
general case of particles with spin s±, s 2 coupling to a total spin S, the angular structure 7 
of a vertex that couples spin (S, m) and orbital angular momentum (£, M — m) to a total 
(J, M) is 

loi = C(S, i, J; m, M - m, M) Y^ m _ M {9, </>) (-l) M ~ m 

with the Clebsch-Gordan coefficients C and a = (J, M, t, S, m, 9, 4>) abbreviates the quan- 
tum numbers and polar and azimuthal angles 9 £ [0, tt] and <p E [0, 2tt]. Note that the 
angles depend on the discretized momenta, 

Tig 

9 = arccos — , <\> = arctan(n J/ /n x ) (2-4) 

\J n l + n l + n l 

with n = (n x ,n y ,n z ) € 1? . For the present case of spinless meson-meson scattering in 
P-wave, 

l* a l a > =3 cos 2 9 (2.5) 

and, thus, eq. (2.2) becomes 

£P-wave = £ 7 * ^ fm) p{q) = £ s j 5 • e ^ ^P-wave = q. (2 g) 

n 

The equation also indicates that for the special case of P-wave, the (5P" wavc function is the 
same as for S'-wave, Gj, given in eq. (2.2) 3 . An analog situation has been found for the 
mass shift of bound states in a box in refs. [82, 83]. This is no longer the case for D- and 
higher waves. Note also that sums ^ n with an S'-wave and a P-wave vertex disappear, 
i.e. there is no mixing of S- and P-waves. As is well known, there is mixing of S- and 
G-waves, and of P- and P-waves, but this is not an issue in the present study because the 
respective higher partial waves are small and, anyway, not generated from the NLO-terms 
used in this study. 

With the G function derived for S- and P-waves, and using the LO potentials V 2 and 
NLO potentials V4, with the low-energy constants of the fit of table 1, eq. (2.3) can now be 

3 This may be seen as follows: in S'-wave, a contribution to the multiplicity at fixed m can be written as 
m = n\ + + n 3- The multiplicity for a set of such rii is 2^ 3_JVz ' 3!/iVg! where N z is the number of zeros 
in (m, ri2, n-j) and N g the number of equal rii. On the other hand, for P-wave, 3 cos 2 9 — in\/{n\ + n\ + 
n 2 ). Weighting each distinguishable permutation of (711,712,71,3) with this factor, it is easy to see that the 
multiplicity at fixed m is the same as for the S-wave case. 
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solved for its tower of zeros and at different box sizes L. The results are shown in figure 3 
for the eight considered reactions. For the (I, L, S) = (0, 0, 0) case, one can perceive the 
avoided level crossing around E = 990 MeV that comes from a combined effect of the KK 
threshold and the /o(980) resonance located there. The situation of resonances close to 
thresholds and methods to extract those resonances from lattice data, e.g. by using twisted 
boundary conditions, have been discussed in refs. [53, 54]. A similar situation is also given 
for the ao (980) [cf. figure 3] although in this case the resonance pole is considerably hidden 
on the pu sheet in the present solution, see table 2. Also in case of the k(800) quantum 
numbers, shown to the upper right in figure 3, avoided level crossing is visible, coming in 
this case solely from the rjK threshold, as the k resonance pole is at lower energies. For the 
P-wave cases of the if* (892) and /?(770) quantum numbers, avoided level crossing is visible 
at higher energies E > 1.1 MeV. This is due to the first free level of the higher channel, i.e. 
7]K for the if* (892) and KK for the p(770) quantum numbers. As these channels couple 
very weakly to the dominant irK and irir channels, respectively, the crossing is almost ideal, 
i.e. the different levels come very close at the crossing energy. Note that, in contrast to 
S-wave channels, there is no avoided level crossing at the respective higher thresholds for 
the P-wave reactions, nor is there a lowest level close to the lower thresholds. An exception 
is given for the ^-channel, but there the level at E ~ 980 MeV arises from the ^(1020) 
resonance that is just below the KK threshold, cf. table 2. The explanation for these 
findings is postponed until section 3.4. 

For the k(800) channel, around the energy of the resonance, there is no structure visible 
in the levels shown in figure 3, and the same applies to the cr(600). For the latter case, 
already studied in ref. [54], it is known that the cr(600) pole can still be reconstructed, 
even if its position far in the complex plane results in a totally washed-out avoided level 
crossing. The situation should be similar for the k(800), and this question is addressed in 
section 3. 

It should be noted that in the present approach we consider the finite volume effects 
on the meson-meson interaction, but do not attempt to determine finite size effects on the 
lattice levels. A comprehensive treatment of such systematic uncertainties would require 
to address the actual lattice action that generates the levels, but that goes far beyond the 
scope of the present study. 

3 Reconstructing the k(800) 

3.1 A suitable potential for resonance extraction 

In the one-channel case, the Liischer formalism provides a model-independent connection 
between phase shift (in the infinite volume limit) and lattice levels [19-21]. For the extrac- 
tion of scattering lengths, see, e.g., also refs. [50, 84]. In the two-channel case (see, e.g., 
ref. [50]), phase shifts, poles and residues can still be reconstructed model-independently 
provided three independent measurements of lattice levels at the same energy, to deter- 
mine the channel transitions Vn, V12, and V^- This information can be obtained, e.g., 
from twisted boundary conditions or using spatially asymmetric boxes, as pointed out in 
refs. [53, 54]. 
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Figure 3. The first few levels for the considered reactions as calculated from the present solution 
for meson-meson scattering. Also, the corresponding resonances are indicated. 
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However, if lattice data are provided with finite accuracy — present-day lattice simula- 
tions do not provide data for more than a few volumes, and with rather large error bars — 
the reconstruction of the infinite volume limit becomes a challenging task. A possible way 
of proceeding is to expand a general scattering potential V in powers of the scattering 
energy, in one or more channels. The expansion parameters are then fitted to the lattice 
levels, and with the same potential V the infinite volume limit can be recovered. Even 
broad resonances like the cj(600) or resonances close to a threshold like the /o(980) can be 
reconstructed as shown in ref. [54]. 

The two-channel fit potential used in ref. [54] is given by 



where i,j = 1,2 indicate the channel. The six parameters Ojj = dji, bij = bji provide 
an expansion up to terms linear in s = E 2 , with the expansion point so = 4Mj^ at the 
KK threshold that is suitable for the study of the /o(980) resonance since the parameter 
correlations are minimized. The problem at hand in ref. [54] is specially suited for the 
choice of the fit potential shown in eq. (3.1), because the LO chiral interaction in the 
(KK, tttt) system depends only linearly on s. If then pseudo-data, generated from the very 
interaction, are fitted with the ansatz of eq. (3.1), ideally a \ 2 °f zero can be obtained. In 
the present case of the (r]K, ttK) system, the LO interaction depends already on s and t, 
and the NLO terms have s 2 -, t 2 ■ ■ ■ -dependence. A possible refinement of the potential V 
could, thus, be an expansion not only in s but also in t, 



The dependence on the Mandelstam variable u can be expressed in this ansatz by s, t and 
a constant, and terms ~ (u — uo) should not be included explicitly. One should note that in 
S-wave, terms proportional to t are to a good approximation linear functions of s, and large 
correlations between a, b, b' appear. It is then more efficient to absorb the t-dependence in 
the parameter b, and also in c for the small order-s 2 differences between s and t. 

As we have just seen, in the presence of terms ~ t and higher order terms from NLO, 
a quadratic term of the form (s — sq) 2 is indispensable in the expansion of the fit potential. 
This will, of course, introduce additional free parameters in the fit potential. Together with 
the fact that only very few data points are available in actual lattice simulations, there 
will still be large parameter correlations and, worse, the fit will have too much freedom in 
between those data points and generate solutions that have a good x 2 but should not be 
accepted because the natural hierarchy of the expansion is no longer given — e.g., the term 
~ (s — so) 2 may be unnaturally large compared to the terms (s — so) 1 and (s — so) . The 
situation can be improved by providing the expansion the well-defined, model- independent 
LO term of the chiral expansion. This term almost saturates the constant and the linear 
parameters a and b, and thus considerably stabilizes the fit. Remaining contributions 
to a and b will then be much smaller, and the above-mentioned, unnatural solutions are 
automatically avoided. 

The question is in which form to cast the resulting fit potential. For a resonance 
analysis, the inverse amplitude method itself is particularly suited because the organization 






(3.2) 
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of the LO and NLO potentials, V2 and V4, allows for solutions in which resonance terms 
can be expressed in terms of the V2, V4 [22]. Thus, this form of the parameterization is 
quite flexible, and we choose it for the final form of the fit potential used throughout this 
study. The finally adopted form of the fit potential reads in the one-channel case 



v = ( T/ 2 ) > Vt = a + b(s - s ) + c(s - s ) z + d(s - s o y + ■■■ (3.3) 



with V2 = Vho the fixed LO term of the chiral expansion. In other words, we take the form 
of the inverse amplitude, leaving the LO term V2 as it is, and expand V4 in powers of s. 
The T matrix in the infinite volume limit is then obtained as 

T= (yzty^G ' Le -> T = V &t + V m GT , (3.4) 

and the lattice levels are given through the condition 

(yfit)-i _ q = o <_> y 2 _ y&t _ y2Q = q ( 3 5 ) 

where G = G n K, G = G n K from eq. (2.2) for the k(800) channel. In the present study, 
we do not generalize this fit potential to two channels. The pseudo-data which will be 
fitted are, of course, generated from the full two-channel solution obtained in the previous 
section. The k(800) is situated well below the r]K threshold, and the influence from the 
r]K channel should be noticed, at most, as a small perturbation that can be absorbed in 
the fit potential. Note in this context that the pseudo-phase [53] in the (KK, tttt) system 
is very close to the true phase up to a few MeV below the KK threshold [54]. In any 
case, these arguments will be tested in the following: if no good fits could be obtained, the 
one-channel extraction method would lead to erroneous phase shifts and pole positions. 

As the LO chiral interaction V2 enters the fit potential explicitly to stabilize the solu- 
tion, the question arises to which extent extracted pole positions are biased: for example, 
V2 alone already generates poles in S-wave that can be identified with the <r(600), /o(980), 
and ao(980) resonances [28]. This is also the case for the k(800) resonance. However, note 
that V^ fit in eq. (3.3) contains also terms ~ s° and ~ s 1 as they appear in V2. Thus, V &t 
is sufficiently general to modify the fixed LO contribution, and resonances may or may 
not appear in the reconstruction process, as dictated by the lattice data themselves. This 
provides a much more model-independent approach to the extraction of the finite volume 
limit as it would be the case by making assumptions on the resonance propagator as, e.g., 
in refs. [4, 9, 15, 17]. 

Finally, it is worthwhile noting that the proposed one-channel extraction method is 
entirely independent of the used regularization scheme because only the difference Gj — Gj 
appears in the connection between the T-matrix in the infinite- volume limit and the lattice 
levels, eqs. (3.4) and (3.5). 

3.2 Extraction of the k(800) 

Fitting pseudo-data with the potential of eq. (3.3) allows to study the propagation of the 
uncertainty from the lattice to the quantities of interest, i.e., the phase shift and the k(800) 
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Figure 4. Left: Pseudo lattice-data and (s°, s , s 2 ) fit to those data with uncertainties (bands). 
Right: Solid line: Phase shift of the solution that was used to generate these pseudo lattice-data. 
Also, the resulting error bands of the (s°, s 1 ), (s°, s 1 , s 2 ), and (s°, s 1 , s 2 , s 3 ) fits to these data are 
shown. 



pole position in the infinite volume limit. To generate the pseudo-data, we consider the 
levels as shown in figure 3 to the upper right. The possible pole position is expected around 
E ~ 650 — 850 MeV, but at the values of L shown in the figure, which are values that can 
be covered in actual lattice simulations, there is no level at these energies. Still, one can 
choose the nearby level at the irK threshold and the following one, without getting too 
close to the rjK threshold. In principle, one could include also data from the r/K threshold 
region, generalizing the fit potential of eq. (3.3) to two channels. However, in actual lattice 
simulations only few data will be available and the number of free parameters in the two- 
channel potential would be three times larger, most of them being weakly constrained (the 
ones from the V^k^k and V^k^k transitions). This would, thus, be a poor strategy to 
study the /c(800). 

Given these arguments, the finally chosen pseudo-data from the levels of figure 3 are 
shown in figure 4 to the left. A constant error of 2 MeV is assigned to each data point, but 
we will also study the case of a 5 MeV error. In ref. [54] the additional error, induced by the 
statistical displacement of the centroids of the error bars, has been studied. Displacing the 
centroid according to a Gauss distribution with half the width of the error bars themselves, 
the uncertainties of extracted phase shifts and pole positions were found to be enlarged by 
~ 40%. We will not repeat this exercise in the analyses of this work, but one should keep 
in mind additional errors of this size on the extracted quantities. 

Using eqs. (3.3) in eq. (3.5), the fit is performed at an expansion point of sq = (750 MeV) 2 
in eq. (3.3). Depending on the power to which the fit potential is expanded, the fits are 
labeled (s°), • • • , (s°, s , s 2 , s 3 ) in the following. Error bands are generated by accepting 
all solutions for which x 2 < Xbest + ^ where Xbest * s the x 2 °f the best fit. The error bands 
for the levels themselves, from the (s°, s , s 2 ) fit, are shown in figure 4 in the left panel, 
and for the phase shift to the right. Also, to the right, the solid line shows the phase shift 
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Figure 5. Pole reconstruction of the k(800). The actual pole position is obtained from a global fit 
to meson- meson scattering data, cf. table 2. Shown are the pole positions of the k with uncertainty 
areas as obtained from different fits to pscudo lattice-data. The lattice data carry an error of 2 MeV, 
except for the case AE = 5 MeV as specified in the figure. 

of the solution that was used to generate the pseudo lattice-data, which lies, as expected, 
in the center of the error bands up to E ~ 900 MeV while it deviates at higher energies. 
This is because, first, no lattice data are fitted beyond E ~ 920 MeV, and second, the 
solution to generate the data is formulated in the full two-channel approach while the fit 
potential is in one channel only. This manifests itself in a cusp of the full solution at the 
T]K threshold, while the extracted phase shifts are smooth there. As the figure shows, fits 
with higher power in the fit potential show a larger spread of possible phase shifts at higher 
energies. This is expected because more free parameters in the fit potential lead to a better 
fit of the level, but also and unavoidably, to a wider error band. 

The next question is to which precision the pole position of the k(800) can be deter- 
mined. The main problem is the large width of this resonance — small uncertainties on 
the real, physical axis at Im E = MeV result in large uncertainties far in the complex 
plane. In figure 5 the pole positions as extracted from the pseudo lattice-data are shown 
for the different fit functions. The central values are indicated with the symbols while the 
areas show the corresponding uncertainties. As for the central values, the s , s 2 ) and 
(s°, s 1 , s 2 , s 3 ) fits are close to the actual pole position, while the (s°) and (s°, s 1 ) fits are 
not. This indicates that the ~ s 2 term in eq. (3.3) is important. As argued in the previous 
section, this term bears part of the t- and u-dependence of the potential, plus parts of 
the NLO interaction. Both turn out to important. Considering the uncertainty areas, we 
observe the same behavior as for the phase shifts: the areas grow with increasing number 
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of fit parameters and higher powers in the potential. In particular, for the (s ) and (s°, s 1 ) 
fits the areas are very small, giving the wrong impression of high precision but the areas 
do not cover the actual pole position. In contrast, the uncertainty area of the (s°, s , s 2 ) 
and (s°, s , s 2 , s 3 ) fits do cover the actual pole position. 

While in our analysis of pseudo lattice-data the actual pole position is known, in the 
analysis of real lattice data it is not. Still, the previous observations lead to the formulation 
of a simple fit strategy: 

1. The fit potential should be expanded to higher and higher powers in s. 

2. Once the central value of the pole position does not change any more — here, given 
for the (s°, s 1 , s 2 , s 3 ) fit — one can assume convergence. 

3. One can then use the previous fit to determine the uncertainty area — here given by 
the s 1 , s 2 ) fit and the yellow area in the plot. This can be regarded as the final 
result of the analysis. 

The chosen data uncertainty of 2 MeV is very small, and we show in figure 5 also the 
uncertainty area for a 5 MeV error, for the (s°, s , s 2 ) fit. This results in an error of 
around ±90 MeV for the real part of the pole position and around ±75 MeV for the 
imaginary part, i.e. ±150 MeV for the width. The residues of the extracted k(800) poles 
and of all other resonances analyzed in this study can be found in section 4.4. 

3.3 Case of a low k(800) pole 

The pseudo-data analyzed in the previous section are based on the global fit to meson- 
meson scattering data carried out in section 2. The outcome for the real part of the k(800) 
pole position of zq = (815 ± i 226) MeV is rather at the upper end of the pole positions 
quoted in the PDG [59]. To cover also low pole positions of the k(800), we repeat the 
analysis using a different solution for which we include, in the x 2 °f the fit, one of the 
lowest pole positions quoted in the PDG [59]. The pole position is the result of the Roy- 
Steiner analysis of ref. [45] with zq = (658 ± i 278.5) MeV. For simplicity, we use only the 
irK channel for the chiral unitary model and include, besides the n pole position, only 
the partial wave data from the ttK partial waves (I,L,S) = (1/2,0,-1), (1/2,1,-1) up 
to E = 1.2 GeV [see figure 1 for the data and their references]. The n pole position is 
included with a chosen error of 35 MeV in the x 2 ■ To constrain the low-energy constants, 
also 7T7T scattering data have been included in the x 2 ■, but only to ensure the <r(600) pole 
is still present. Furthermore, we have made sure that the fitted combination of low-energy 
constants leads to qualitatively correct solutions in the other meson- meson reactions. 

The refit of the low-energy constants (except for Lj that does not contribute in 
irK — > irK) leads to a quantitative description of the data in the irK partial waves 
(I, L, S) = (1/2, 0, —1), (1/2, 1, —1). This is shown by the solid line in figure 6 in the right 
panel (the K* (892) channel is not shown again). The pole position at zq = (679± i 271) MeV 
is quite close to the fitted one of ref. [45] and at the lower end of values quoted in the 
PDG [59]. 
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Figure 6. Case of a low k(800) pole. Left: Pseudo lattice-data and uncertainty of the (s°, s 1 , s 2 ) 
fit to those data (bands). Right: Solid line: Phase shift for the solution that was used to generate 
these pseudo lattice-data. Also, the resulting error band of the (s°, s 1 , s 2 ) fit to these data is shown. 
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Figure 7. Case of a low re(800) pole. Pole positions of the k (central values) together with 
uncertainty areas, from fits to pseudo lattice-data. All symbols and areas are defined as in figure 5. 

The fit solution provides the pseudo-data that are selected following the same argu- 
ments as in the previous section. The pseudo-data and the error bands of the corresponding 
(s°, s 1 , s 2 ) fit are shown in figure 6 to the left, applying the same Xbest + 1-criterion as in 
the previous section. To the right, the error band of the extracted phase shift from the 
(s°, s 1 , s 2 ) fit is indicated with the filled (cyan) area. Comparing with the previous case 
of a higher K-pole, shown in figure 4, the error bands both for levels and phase shifts are 
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Figure 8. Levels 1 and 2. Thick (red) solid lines: case of higher pole position from section 3.2. 
Thin (black) solid lines: case of lower pole position from section 3.3. Dashed line: Approximation 
of cq. (3.7) for the case of lower pole position. The squares to the left indicate the measurements 
taken to calculate the approximations to phase shifts and pole positions shown in figure 9. 

similar for the respective (s°, s , s 2 ) fits. 

In figure 7 the pole positions as extracted from the pseudo lattice-data are shown for 
the different fit functions, in full analogy to the results shown in figure 5. The central 
values are indicated with the symbols, the areas show the corresponding uncertainties. We 
can perceive the same behavior of the fits as in the previous case of a rather high K-pole 
position: expanding the fit potential to too low powers, wrong central values are extracted 
in combination with too small uncertainty areas that give the impression of an unrealistic 
precision. Once also the s 2 -term is fitted, the central value is close to the actual pole 
position. Including also the s 3 -term in eq. (3.3) serves then to confirm the result of the 
(s°, s , s 2 ) fit, because the central value is not changed any more and only the uncertainty 
area grows. The realistic uncertainty is then again given by the yellow shaded area from 
the (s°,s\s 2 ) fit. 

3.4 Qualitative analysis using the lowest level 

In the previous sections we have seen that relatively precise lattice data are needed to pin 
down the pole position of the k(800). Still, one may wonder if with a few and low precision 
data points from the lowest level, which is the one that can be most precisely determined in 
actual lattice simulations, one can at least make qualitative statements about the presence, 
absence, and maybe approximate position of the pole. 

In figure 8 the levels 1 and 2 are shown. The case of a higher k(800) pole position [cf. 
section 3.2] is indicated with the thick (red) solid line, the case of lower pole position [cf. 
section 3.3] with the thin (black) solid line. As the figure shows, the largest differences of 
the two solutions are visible for level 1, and those differences are enhanced to more than 
10 MeV for the smallest L. In contrast, level 2 shows only differences of a few MeV below 
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E ~ 900 MeV. In the following, we will study to which extent the low-L enhancement of 
level 1 can help to draw qualitative conclusions on the k(800) pole. 

We start with the derivation of an approximate expression for the L-dependence of 
the first level. In the present formulation the function G of eq. (2.2) can be expanded in a 
Laurent series around the pole at threshold, 



G 



1 



Ami m 2 L 3 E — (mi + m 2) 



+ 0[E - (mi + m 2 )] c 



(3.6) 



where we can ignore contributions from G and Gs in eq. (2.2) of order 0[E — (mi+ tt^)] 1 / 2 . 
For an S- wave interaction, we can approximate the potential V close to threshold by a 
constant, V &t = C^k- The levels are given by 1 - V &t G = (the problem can be reduced 
to one channel for this exercise), so that 



E = mi + 7712 + 



K 



Ami m 2 



(3.7) 



This is equivalent to the standard result for the scattering length, see refs. [19, 20, 84]. 
There will always be a level close to threshold with this functional behavior, because G 
diverges at threshold and always equals (V^)" 1 for some E close to it. In particular, 
we read from eq. (3.7) that if the interaction is attractive (C n K < 0), the level will be 
slightly below threshold, while for C n x > it will be above. For the 5-wave, we observe 
exactly this behavior for all cases shown in figure 3. Note that once C n x is determined 
from lattice data using eq. (3.7), one can relate it to the scattering length. With the 
isospin I = 1/2 5"-wave scattering length defined as lim gc m _,.o <?c.m. cot<^y 2 = V a i/2 an< ^ 
-ReT® /2 (E = M n + M K )/[8ir(M n + M K )\ in the current normalization of the 



6 l/2 



amplitude T we obtain 



1 



A 



h/2 



87T E C W K G-kK — 1 



(3.8) 



E=hU+M K 



Back to the original problem, the approximation of eq. (3.7) for level 1 is shown in 
figure 8 with the dashed line for the case of low k pole. The agreement is quite good. With 
C n x determined from eq. (3.7), we can exploit the fit potential of eq. (3.3). The measured 
C n K fixes the parameter a of V^ 1 , while b, c, d are set to zero as they cannot be fixed, 



T/thres. 
v 2 



ythres. 



const., 



ythres. 



V 2 {b = (M v + M k ? 



(3.9) 



Using eq. (3.4), phase shifts and pole positions may be determined now. For the numerical 
calculation, we have determined C^k by measuring the lowest level at L = 2.5M~ 1 , as 
shown by the squares in figure 8, and then using eq. (3.7). The result is shown in figure 9 
for the two cases of higher and lower pole position. In the left panel, the actual phase shifts 
(solid lines) and the ones determined using eq. (3.9) are shown (dashed lines). Indeed, using 
only the information from the lowest level, encoded in C^k, the phase can be qualitatively 
described, in particular for the case of a higher k pole. The pole positions are shown in the 
right panel of figure 9. The actual pole positions are indicated with the squares, the ones 
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Figure 9. Dashed lines/dashed circles: Approximations of phase shifts (left) and pole positions 
(right) using only one data point from level 1 as shown by the squares in figure 8. For comparison, 
the full results are shown: Thick solid (red) line and thick (red) square: case of high k pole position. 
Thin solid (black) line and thin (black) square: case of low n pole position. 

determined using eq. (3.9) are shown with the circles, and the pole position obtained when 
setting V^ fit = 0, i.e. the pole position from the LO calculation (in one channel), is marked 
with the cross. Here, the situation is not as clear as for the phases. In particular, the 
approximations of the actual pole positions, given by the method proposed here (circles), 
show no systematic improvement compared to the case V^ fit = (cross). In summary, 
measuring F 4 fit at one L using the lowest level improves the description of the phase close 
to threshold. However, beyond the threshold region, clearly the higher order terms ~ s, s 2 
in V^ fit need to be known to allow for a quantitative determination of the phase and the k 
pole. 

To conclude this section, we come back to a question left open in the discussion of 
figure 3: why is there no level close to the lower thresholds of the P-wave reactions, while 
there is always one for the S- wave reactions? For P- waves, the potential close to threshold 
will be of the form V = a' f/c. m . with a nearly constant a'. Expanding q c . m . as a function of 
E one finds that the pole of G (see eq. (3.6)) is precisely canceled in 1 — VG = 0, 

1 - VG ~ 1 - - r „ . a r (3.10) 

2L 3 (mi+m 2 ) 

and there is no need for this expression to become zero, i.e. no need for the occurrence of 
a level close to a P-wave threshold. The absence of avoided level crossing at the respective 
upper thresholds, for the P-wave channels, can be explained similarly. 

4 Extraction of other excited meson states 

To put the previous results on the k(800) into context, in particular the required high pre- 
cision of the lattice data, in the following the cr(600) is analyzed that has similar properties. 
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Figure 10. Left: Pole reconstruction of the <r(600). The (black) symbols show different fits, 
the actual pole position is shown by the (red) open square. Uncertainty areas are labeled. For 
comparison, the result from ref. [54] is shown with the hatched circle (pole position) and hatched 
uncertainty area, obtained from a (s°, s 1 ) fit. Right: 6 data points from the upper level (level 2) 
plus 5 points from level 1, as used for the fits. The data errors and uncertainty band for one fit is 
shown as indicated in the figure. 

Also, the results for these S"-wave resonances are contrasted with the P-wave K* (892) and 
p(770) states, as well as the repulsive irK and 7T7T phase shifts. 

4.1 The er(600) at next-to-leading order 

In the previous sections, lattice data with an error of 2 MeV were required to obtain 
a sensible error on the pole position of the k(800), as shown with the (yellow shaded) 
(s°, s , s 2 )-areas in figures 5 and 7. This result is in certain contrast with the pole recon- 
struction of the cr(600) in ref. [54] where a few data points with a 10 MeV error were found 
enough to obtain a reasonable determination of the pole. We repeat the current analysis 
to put these differences into context. For this, we use the same values of L as in ref. [54] 
to generate 6 data points with a 10 MeV error 4 . These data are fitted in a (s , s 1 ) fit, 
like in ref. [54], with the result for the pole position and uncertainty (orange area) shown 
in figure 10 to the left. The corresponding result from ref. [54] is shown with the hatched 
areas in figure 10. As the figure shows, for the present result and the one from ref. [54], 
the uncertainty is of quite similar size. However, the central value in the present case 
(filled circle) is quite far from the actual position: Like in the case of the ac(800), there 

4 To generate the data, the global fit from section 2 needed to be slightly modified, because that solution 
contains, in the (I, L, S) = (0, 0, 0) quantum numbers, a pole in the T-matrix below the nir threshold at 
E — 200 MeV. Note that also the solution of ref. [22] has a pole below threshold for these quantum numbers. 
This is not noticeable in the phase shift, but the first lattice level, below the irn threshold, is affected at 
small L. Such effects close to threshold are, apart from being unphysical, difficult to accommodate in the 
expansion of the fit potential. 
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are now sizeable NLO terms ~ s 2 that cannot be accommodated in a simple (s°, s 1 ) fit. 
Thus, to better reconstruct the original pole position of the <r(600), we have performed an 
(s°, s , s 2 ) fit, but with a 2 MeV error for the data (yellow area). As the figure shows, the 
corresponding uncertainty area (yellow) is of similar size as for the case of the k(800) from 
Figs. 5 and 7. Concluding, in ref. [54] a 10 MeV error on the data produces an uncertainty 
of similar size as 2 MeV data errors here, because the fit potential here is expanded to 
higher powers. As figure 10 shows, even the (s°, s , s 2 ) fit shows no agreement of the cen- 
tral value with the actual pole position. The (s°, s , s 2 , s 3 ) fit, however, does. To minimize 
the uncertainty area for this case, in addition to the 6 points from level 2, 5 points from 
level 1 are included in the fit, that would be known in an actual lattice calculation anyway. 
As the figure shows, one can then provide again a 10 MeV error to the data to obtain a 
large, but still limited uncertainty area 5 . 

4.2 P-wave states: Extraction of the K*(892) and the p(770) 

While the broad k(800) and o"(600) S-wave resonances, discussed in the previous sections, 
require relatively precise lattice data, some P-wave resonances are expected to show clearer 
signals on the lattice: For P-wave the lowest lattice level — i.e. the most accessible one — 
is located already in the resonance region, in contrast to the S-wave that has a level at 
the lowest threshold (see figure 3 and section 3.4). Also, resonances like the K*(892) and 
p(770) are much narrower than the analyzed 5-wave states and the corresponding ttK and 
7T7T channels show only very weak coupling to the respective r]K and KK channels [22], 
making the analysis easier than in case of, e.g., the /o(980) or ao(980) resonances [53, 54, 58]. 
Lattice calculations and resonance extraction at finite volume for the p(770) have already 
allowed for a determination of its mass and the pnir coupling constant, at a pion mass 
of around twice the physical one, see, e.g., ref. [15]. The relativistic Breit-Wigner form, 
assumed in ref. [15], or effective range formulae assumed in refs. [4, 9], are adequate for the 
p(770), but the method presented here is a more general model-independent extraction of 
the resonance properties. 

The S-wave fit potential from eq. (3.3) can be generalized to P-wave with minor 
changes. Noting that the threshold behavior is given by Vpl wave ~ <2c.m.> we can simply 
multiply this factor to the S'-wave potential to obtain 



^P-wave = ' / , Vf = <Zc.m. [« + &(* " «o) + c(s - S f + ■ ■ ■ ] (4.1) 



and V% is given by the P-wave projected LO term. 

For the generation of pseudo-data in the (r/K, ttK) P-wave channel, the global solution 
determined in section 2 has been used, while for the pseudo-data of the (KK, tttt) P-wave 
we use the solution of ref. [22] that is slightly better for this quantum number (dashed vs. 
solid line in figure 1). 

5 Note that according to the fit strategy formulated in section 3.2, in an analysis of lattice data where 
the actual pole position is unknown, one would still have to carry out the (s°, a 1 , s 2 , s 3 , s 4 ) fit to confirm 
the pole position found in the (s°, a 1 , s 2 , s 3 ) fit and accept the corresponding uncertainty area as the final 
result. As we know the actual pole position for the pseudo-data, there is no need to perform this task here. 
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Figure 11. Extraction of the if* (892) resonance using 7 data points from level 1 with 10 MeV 
error (left). Center and right: Phase shifts and pole position, respectively. The (yellow) areas show 
the uncertainty of the (s^s 1 ) fit, the (gray) areas the one of the (s°, s 1 , s 2 ) fit. The actual phase 
[pole position] is indicated with the thick solid line [open (red) square] . The central values of the 
fitted pole position are indicated with the filled circle and the filled square. 
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Figure 12. Extraction of the /?(770) resonance using 5 data points from level 1 with 10 MeV error 
(left). Center and right: Phase shifts and pole position, respectively. Labeling as in figure 11. 

The analysis of the pseudo-data can be carried out as in the previous sections, with 
the result for the if* (892) shown in figure 11 and for the p(700) in figure 12. For both 
resonances, the (s°,s 1 ) fits lead to pole positions very close to the actual ones, indicated 
with the open (red) squares. Note that contributions to the s 1 -term in P-wave correspond, 
among others, to terms of the form ts, us, t 2 , or u 2 as they appear in the NLO potential [22] . 
The (s°, a , s 2 ) fits return practically the same pole positions and much enlarged uncertainty 
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Figure 13. Extraction of the repulsive irK (left) and irir phase shifts (right). The thick solid lines 
show the actual phase shifts. The (yellow) bands show the extracted phase shifts using one data 
point of the respective lowest level, at L = 2.5 M" 1 and with a 10 MeV error. 

areas. This indicates that the (s ,^ 1 ) fits and their uncertainty areas can be regarded as 
the final results, according to the fit strategy formulated in section 3.2. 

The observed convergence for (s°, s 1 ) fits should be seen in comparison to the previous 
results for the S-wave resonances: for the /c(800), convergence occurs only for the (s°, s , s 2 ) 
fit, while for the cr(600) the potential needs to be expanded to even one more power. Thus, 
the properties of the considered P-wave resonances are indeed much easier to determine 
than their S"-wave counterparts. 

For the if* (892), seven data points up to values of L = 4.5 M~ l have been fitted. As 
figure 11 shows, the (yellow) uncertainty area of the levels is broad at low L and narrow 
at large L. This indicates that the data at small L predominantly fix the fit result, while 
those at high L may be dispensable. Thus, for the analysis of the p(770) the two data 
points at L > 3.5 M" 1 were omitted as shown in figure 12. Indeed, the p(770) properties 
could be determined to a similar quality as the ones of the if* (892), taking into account 
that the larger width of the p(770) automatically leads to an enlarged uncertainty area. 

Thus, the analysis of pseudo-lattice data can help to identify regions of the lattice spec- 
trum that are particularly sensitive to resonance properties. Such findings can contribute 
to minimize the required effort of actual lattice calculations. 

4.3 The ttK and irir repulsive phase shifts 

To put the resonance analyses of the previous sections into context, we apply the same 
method to phase shifts without resonances. The repulsive ttK isospin 3/2 and 7T7T isospin 
2 phases show no structure and an almost linear energy dependence. As the tttt isospin 
2 phase shift provides a channel with maximal isospin, it is more accessible in lattice 
calculations [10, 13] because disconnected diagrams are absent (the same holds for the 
isospin 3/2 irK S-wave). As figure 3 shows the second level is already at quite high 
energies, for realistic values of L and physical pion masses. Here, we restrict the analysis 
to the lowest level, i.e. we proceed as in section 3.4. One could use several data points from 
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Figure 14. Residues of the resonances. The labeling of symbols and areas is precisely as in the 
figures of the respective resonances, 5,7, 10, 11, and 12. In the left panel, both cases of a heavy 
k(800) from section 3.2 and a light k(800) from section 3.3 are shown. The central panel shows the 
residues of the cr(600) discussed in section 4.1, and in the right panel the residues of the if* (892) 
and p(770) from section 4.2 arc displayed. 



the lowest level, but, as shown in eq. (3.7), the L-dependence is quite fixed and one data 
point at low values of L is enough. For the analysis, the data point is taken at L = 2.5 M" 1 
with an assigned error of 10 MeV, using the values for the low energy constants of ref. [22] 
that produce a slightly better solution for the repulsive phases than the global one derived 
in this study (dashed vs. solid lines in figure 1). As discussed in section 3.4, using one data 
point allows only for a (s°) fit of the fit potential of eq. (3.3). The results for both repulsive 
phase shifts are shown in figure 13. The actual phase shift lies inside the uncertainty area 
up to energies of 500 MeV above the ttK threshold and 400 MeV above the ittt threshold, 
respectively. Thus, due to the almost linear energy dependence of the phases, a single 
data point is enough to reconstruct them up to quite high energies above the respective 
thresholds. 

4.4 Residues 

To finish the resonance analysis, the residues of all analyzed resonances with their uncer- 
tainties are shown in figure 14. The displayed residues are defined as 

T = + 0(E - z )° (4.2) 

& — zq 

where zq is the pole position in the upper half plane and T is the transition within the 
lighter channel (on the second Riemann sheet), i.e. ttK — > ttK for the k(800) and the 
if* (892); vrvr -> vrvr for the <t(600) and the p(770). Note that the residues for the poles 
at Zq in the lower half plane are given as o_i(zq) = al 1 (zo)- For the values of the actual 
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residues, shown by the open (red) squares, the amplitude T in eq. (4.2) is given by eq. (2.1), 
i.e. the solution that was used to generate the pseudo-data. For the fitted residues (central 
values, shown by filled (black) symbols, and areas), T is given by eq. (3.4), using eq. (3.3) 
or (4.1) for the fit potential, depending on whether S- or P-wave states are analyzed. 

The results shown in figure 14 suggest a similar interpretation as for the respective pole 
positions shown in figures 5,7, 10, 11, and 12: The residues of the broad 5-wave resonances 
are much more difficult to determine than the residues of the P-wave states and require 
higher precision lattice data. For example, a 5 MeV error for the lattice data of the k(800) 
barely fixes the sign of the real part of the residue. In contrast, with only a few lattice 
points from the lowest level and 10 MeV error, the (s°, s 1 ) fits of the K*(892) and />(770) 
resonances lead to a very narrow error band on the respective residues. 

5 Conclusions 

For the reconstruction of the k(800) resonance from lattice data, we have first performed 
a global fit to meson-meson partial wave data using the inverse amplitude method based 
on next-to-leading order terms of the chiral expansion. This allows to determine the low 
energy constants and to predict a pole position of zq = (815 + i 226) MeV for the k(800). 
This solution was used to generate pseudo-lattice data that were subsequently analyzed 
expanding a suitable fit potential systematically in powers of s. The key point is that 
the lowest-order chiral interaction serves as explicit input to the fit potential, considerably 
stabilizing the results. With a simple fit strategy the ttK phase and k(800) pole position 
together with the respective uncertainties could be determined in a model-independent 
way. To cover the possible pole positions of the k(800), we have also analyzed the case of 
a lighter k(800). As has been shown, a comprehensive resonance analysis requires rather 
high-precision lattice data. Still, using only one data point from the lowest level, one can 
at least make qualitative statements about the k(800). To put these results into context, 
we have also analyzed other channels using the same technique, showing that the P-wave 
resonances K*(892) and p(770) are much easier to extract than the S'-wave k(800) and 
cj(600) resonances. The proposed analysis method allows not only to quantify the required 
precision of lattice data, but also helps to identify regions of the lattice spectrum that 
are particularly sensitive to resonance properties, which can contribute to minimize the 
practical effort in actual lattice simulations. 
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